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We report results of inelastic neutron scattering measurements of phonon dispersions in optimally 
doped YBa2Cu3 06.95 and compare them with model calculations. The focus is on the in-plane 
oxygen bond-stretching phonon branches. The study of these modes is complicated by anticrossings 
with c-axis-polarized branches; such effects are interpreted through lattice-dynamical shell-model 
calculations. The in-plane anisotropy of the bond-stretching phonons was firmly ascertained from 
measurements on a detwinned sample. Studying not only the in-plane modes involving in-phase 
motion for the two Cu-O layers within a unit cell but also those with opposite-phase motion was 
of great help for establishing a clear experimental picture. The measurements confirm that the 
in-plane oxygen bond-stretching phonon branches disperse steeply downwards from the zone center 
in both the a and the b directions indicating a strong electron-phonon coupling. For the b-axis- 
polarized bond-stretching phonons, there is an additional feature of considerable interest: a sharp 
local frequency minimum was found to develop on cooling from room temperature to T = 10 K at 
wave vector q~0.27 r.l.u.. 

PACS numbers: 74.25.Kc,63.20.Kr,74.72.Bk 
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I. INTRODUCTION 

Numerous inelastic neutron scattering investigations ^ 
0, El on high-T c superconductors and their insulating 
parent compounds have shown that the high-energy Cu- 
O bond-stretching modes soften considerably on dop- 
ing. It was advocated from the very beginning^] that 
this phonon softening indicates a strong electron-phonon 
coupling, but it attracted nevertheless only limited at- 
tention because phonons were widely considered to be 
irrelevant for the mechanism of high-T c superconduc- 
tivity. Recently, renewed interest in the phonons came 
from observations made by angular resolved photoemis- 
sion spectroscopy of an abrupt change of electron velocity 
at 50-80 meV which was interpreted as evidence for an 
ubiquitous strong electron-phonon coupling in high-T c 
superconductors [jj, |(| . Further interest was generated by 
theories!]], H,0| predicting an inhomogeneous charge dis- 
tribution in the Cu-0 planes, in particular in the form of 
stripes. In superconducting samples, charge stripes are 
assumed to be dynamic in nature which makes it very dif- 
ficult to detect them. They m ight reveal themselves by 
coupling to the phonons jlOt [III Il2fl thus making inelastic 
neutron scattering an appealing technique. Naively, one 
might expect precursor phenomena to charge-density- 
wave (CDW) formation similar to what has been ob- 
served in one-dimensional conductors such as KCP some 
time ago|l3j. There, a dip develops in the phonon dis- 



persion of the longitudinal acoustic branch in a narrow 
range of momentum transfer about a wavevector related 
to the spanning vector 2kp of the Fermi surface. 

In the cuprates, the most pronounced effects are ex- 
pected for the Cu-0 bond stretching vibrations because 
the Cu-0 bond length will be significantly modulated 
by the extra charge residing on the stripes. Early 
investigations |14j of the phonons in YBa2Cu307_ :c have 
indeed shown that the O bond-stretching vibrations 
within the Cu-0 planes behave in a normal way in in- 
sulating samples but in an apparently anomalous way in 
optimally doped ones: the frequency of these phonons de- 
creases abruptly when going from the zone center along 
the [100] (=a*) direction or the [010] (=b*) direction 
about half way to the zone boundary. In addition, the 
phonon peaks become very poorly defined in this q-range. 
Unfortunately, twinning of the samples did not allow one 
to decide whether the phonons in question really acquire 
a very large intrinsic linewidth or whether the observed 
peaks are so broad simply because the superimposed con- 
tributions from the a* and the b* directions differ in en- 
ergy. A recent inelastic neutron scattering investigation 
|T5| using pulsed neutrons came to the conclusion that 
the Cu-0 bond-stretching branch along the b* direction 
has a steep, but continuous downward dispersion whereas 
the branch along a* splits into two in the middle of the 
zone. The latter finding was tentatively interpreted as a 
signature of charge stripes. We will discuss this interpre- 
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tation in the light of our own results. 

The present paper describes results of inelastic neu- 
tron scattering measurements made by the triple-axis 
technique on detwinned as well as on twinned sam- 
ples of optimally doped YBa2Cu3C>6.95 aiming at elu- 
cidating the anomalous behavior of the bond-stretching 
modes. The measurements on the detwinned samples 
were key to achieving a better understanding of the 
a-b anisotropy and of the seemingly anomalous line- 
shapes. We show that very complex intensity distribu- 
tions can be explained as resulting from an interaction 
of phonon branches having the same symmetry but a 
different polarization^^. The line broadenings due to 
anti-crossings partly mask a more important reason for 
anomalous lineshapes, i.e. a sharp drop of phonon fre- 
quencies within a narrow range of wave vectors along b* 
observed at low temperatures. As has been pointed out 
in a separate publication|l7| this peculiar behavior is ob- 
served only at low temperatures and strongly suggests 
dynamic charge stripe formation with a period of about 
four lattice constants. 

The rest of the paper is organized as follows. Exper- 
imental details are given in the following section. The 
model used to simulate the phonon dispersions and scat- 
tering intensities is described in Sec. III. These calcula- 
tions are essential for making sense of the measurements, 
which are presented in Sec. IV. A discussion of the re- 
sults and conclusions are given in Sec.V and VI, respec- 
tively. Appendix A provides a brief description of anti- 
crossing behavior between neighboring phonon branches. 
The special force constants, beyond the shell model, re- 
quired to simulate the observations are described in Ap- 
pendix B. 



II. EXPERIMENTAL 

The first series of measurements was performed on 
a detwinned sample. It consisted of 26 detwinned sin- 
gle crystals co-aligned with an effective mosaic spread 
of ~3 degrees. The total volume of the sample was ~1 
cm 3 . Unfortunately, the crystals cleaved into many small 
pieces after the first cooling cycle, and therefore addi- 
tional measurements had to be performed on twinned 
samples. The next round of measurements was carried 
out on a composite sample weighing a total of 31 g which 
was used previously for a par t of the pulsed neutron mea- 
surements of Chung et al|l5|. However, this sample was 
not very well suited for triple-axis measurements because 
the lateral dimensions were too large compared to the 
size of the neutron beam. Moreover, high-precision mea- 
surements of the c-axis lattice constant of the individual 
crystals revealed that some of them were underdoped. 
Based on these findings, another composite sample was 
assembled out of 3 optimally doped (YBa2Cu3 06.95, T c 
= 93 K) single crystals of combined volume of ~1.5cm 3 
co-aligned with a mosaic spread of 2.2 degrees. All the 
temperature dependent studies were done on this sample. 



The experiments were carried out on the triple-axis 
spectrometer IT located at the Orphee reactor using dou- 
bly focusing monochromator (Culll or Cu220) and an- 
alyzer (PG002) crystals. Specifically, Culll was used as 
monochromator for the measurements on the detwinned 
sample to maximize the intensity whereas Cu220 was 
used for the measurements on the twinned samples to 
achieve high resolution. The energy resolution at energy 
transfer of E«70 meV was AE=4.6 (2.9) meV (FWHM) 
for the Culll (Cu220) monochromator and a final en- 
ergy E^=14.7 meV. The measurements were carried out 
in different scattering planes in order to elucidate the 
eigenvectors of the atomic vibrations in question. Most 
measurements were performed at T=12 K but selected 
phonons were also studied at higher temperatures up to 
T=300 K. 



III. PHONON MODEL CALCULATIONS 

In measurements on compounds with many atoms in 
the unit cell, an assignment of phonon peaks to partic- 
ular modes has to be based on model calculations. For 
this purpose, extensive model calculations were carried 
out prior to and in parallel with the experiments. The 
model used as a starting point was the common interac- 
tion potential model reported in Ref. 19, which is quite 
successful in describing the phonon dispersion curves of 
a number of cuprates. In the framework of this model, 
YBa2Cu307 is considered as an ionic compound, and 
the interatomic interactions are modeled as a sum of 
Coulomb forces and short range repulsive forces. In addi- 
tion, the polarizability of the atoms is taken into account 
using the shell model formalism. For a metallic com- 
pound like YBa2Cu307 a term accounting for screening 
by free carriers is added. This model gives a decent de- 
scription of most of the available data shown in Ref. 1. 
We further improved upon the agreement between model 
and experiment by tuning the parameters; the refined 
parameters are listed in Table I. However, it was already 
clear before the present study that extra ad hoc inter- 
actions are required to adequately describe the oxygen 
bond-strctchingmodes . 

Experiments^], 0] have shown that the bond- 
stretching phonon frequencies are strongly renormalized. 
This effect is much stronger in the [100] and in the [010] 
directions than in the [110] direction.lt leads to a down- 
ward dispersion from zone center to zone boundarv|l4j 
which cannot be reproduced by simply tuning the pa- 
rameters of the common interaction potential model[18j. 
Adding to the model a negative breathing deformabil- 
ity lowers the calculated energy of the planar breath- 
ing mode, which is the endpoint of the bond-stretching 
phonon branch in the [110] direction. At the same time, 
this term lowers the energy of the corresponding branches 
in the [100] and the [010] directions, but by only half the 
amount of that in the [110] direction. Therefore, addi- 
tional terms were added lowering specifically the energies 
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of the " half- breathing" modes at the zone boundary of 
the [100] and - with a different parameter - of the [010] 
direction. The extra terms are described in detail in Ap- 
pendix B. 

To proceed further with this discussion, it is neces- 
sary to illustrate the situation with actual data. Fig- 
ure 1 compares calculated dispersions with measured O 
bond-stretching modes dispersing along [010] and involv- 
ing opposite-phase motion for the two Cu-0 layers within 
a unit cell (A4' symmetry). The fit to the 200-K data 
with the additional terms discussed above is indicated 
by the magenta curve. This fit matches the data well 
at the zone center and boundary, but falls significantly 
above the data in the middle of the zone. To remedy this 
deficiency, a further term was included that produces a 
maximum effect at the halfway point but zero effect both 
at the zone center and at the zone boundary. The result- 
ing fit is indicated by the red line in Fig. 1. An impor- 
tant result of the present investigation is the observation 
of an anomalous softening in the [010] direction within 
a rather narrow range of wave vectors on cooling from 
room temperature to low temperatures. In order to fit 
the 10 K data, an additional term producing a sharper 
effect halfway to the zone boundary had to be included 
with the resulting fit indicated by the blue line in Fig. 1. 
We believe that all of the extra terms we have added re- 
flect the electron-phonon coupling resulting from doping 
holes into the CuC>2 planes. 

Before continuing on, it may be useful to make 
contact with recent theoretical work. Bohnen, Heid, 
and Krauss |2^] have calculated the phonon dispersions 
for YBa2Cu3C>7 using a method based on evaluating 
the band structure in the local-density approximation 
(LDA). Their results successfully match the experimen- 
tal A 4 ' branch at 200 K (see Ref. [13) though they do 
not predict the temperature induced softening halfway to 
the zone boundary. Further, their calculations reproduce 
the corresponding A4 branch along [100] very well (Fig. 
8). The LDA results provide some justification for the 
ad hoc interactions that have been discussed above and 
point at the unconventional nature of the temperature 
dependence observed in the [010] direction. 

We note again that Fig. 1, which illustrates the effects 
of the special terms, was prepared for branches of A4- 
symmetry. Phonons of A4-symmetry are distinguished 
from those of Ai-symmetry in that the atoms in the Cu- 
O bi-layer move in opposite phase instead of in-phasc. 
For modes of Ai-symmetry, the downward dispersion 
of the bond-stretching mode leads to anti-crossings with 
branches of the same symmetry which originally have a 
different polarization but that hybridize with the bond- 
stretching phonons where they are close in energy (read- 
ers who are not familiar with such effects are referred 
to APPENDIX A). The anti-crossings are naturally ex- 
plained by the model, and the model calculations for Ai 
frequencies and intensities agree well with the experimen- 
tal measurements, as will be shown in the next section. 
It happens that the branches of A4-symmetry are much 



less affected by anti-crossing effects. 



IV. RESULTS AND ANALYSIS 
A. Low temperature measurements 

Experimental spectra obtained on the detwinned sam- 
ple are plotted in Figs. 2,3. The two series of scans were 
conducted along lines from the (4,1,0) reciprocal lattice 
point to (3.5,1,0) and from (1,4,0) to (1,3.5,0), respec- 
tively. These lines were chosen because (i) the inelastic 
structure factors of the bond-stretching modes are quite 
favorable here and (ii) these choices allowed us to ex- 
ploit focusing effects. The figures further show simulated 
spectra calculated with the help of the lattice dynamical 
model described in the previous section. We note that the 
full resolution in q and w was taken into account for the 
simulations. There is obviously very good agreement be- 
tween the experimental and the calculated spectra. An- 
other way to visualize this agreement is to represent both 
types of spectra as (color coded) contour plots (Figs. 
4,5). The plots computed from the simulated spectra 
also show the phonon dispersion curves calculated from 
the same model for the high symmetry directions. All 
the depicted branches are of the same symmetry, i.e. Ai 
(Figs. 2,4) and Ai'(Figs. 3,5), respectively (in order to 
distinguish phonon branches along a* from those along 
b* we add a prime to the symmetry labels for the b-axis 
branches). They are sometimes simply labelled as longi- 
tudinal optic branches. However, this labelling is rather 
imprecise as there are other phonon branches with a dif- 
ferent symmetry having a longitudinal character as well. 
On the other hand, the polarization may be very differ- 
ent for different branches of the same symmetry class: 
as will be discussed below, some of the Ai(Ai')branches 
are c-axis polarized and hence are not observable for mo- 
mentum transfers within the [100]-[010]-scattering plane. 
Nevertheless, they cannot be ignored in the interpreta- 
tion of our data because they hybridize with in-plane 
polarized branches in certain parts of the Brillouin zone 
(because they belong to the same symmetry class) and 
hence show up in the scans under discussion. 

Our conclusions regarding c-axis-polarized Ai 
branches are based on previous studies. In particular, 
the c-axis polarized branches with zone-center energies 
of 62 meV and 54 meV, respectively, were studied 
in Ref. 14, where they were measured with a large 
c*-component of the momentum transfer. A summary 
of the available data from the present and from previous 
studiesflj together with the calculated dispersion curves 
in the whole range of phonon energies is shown in Fig. 
6. 

As already indicated in Sec. II, a very important 
step in establishing the dispersions of the bond-stretching 
modes has come from studies of branches with A4- 
symmetry, i.e. the modes where the atoms in the Cu- 
O bi-layer move in opposite phase. We found that the 
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bond-stretching vibrations of Ai-symmetry and those of 
A^symmetry have very similar energies, as one should 
expect from the weak coupling across the Y atom. What 
makes a study of the A4-bond-stretching branches very 
attractive is the fact that mixing with c-axis polarized 
modes is less of a problem. The A4-modes had to be 
studied on a twinned sample. The measurements were 
made around the (3,0,2) reciprocal lattice point. Repre- 
sentative spectra are shown in Fig. 7; see Ref. 17 for 
further data. As one can see, high resolution allowed us 
to separate the modes for the a* direction and for the 
b* direction, respectively. The results for the a* direc- 
tion are summarized in Fig. 8. The results for the b* 
direction are included in Fig. 1. 

The good agreement between our model and the exper- 
imental gives support to a rather conventional picture of 
phonon dispersion with a notable exception of anoma- 
lous dispersion of plane oxygen bond-stretching vibra- 
tions calling for the special terms. The following pic- 
ture emerges from our analysis: the branches starting at 
E=67 meV and 73 meV, respectively, have in-plane Cu-0 
bond-stretching character with the atoms in the Cu-0 bi- 
layer moving in-phase. The energy difference at the zone 
center can be explained entirely by the difference in the 
Cu-0 bond lengths along a and b. Both branches show a 
steep downward dispersion half-way to the zone bound- 
ary. The b-polarized branch even has a pronounced local 
minimum at q«0.27 r.l.u. whereas the a-polarized branch 
does not have such a local minimum but is quite flat for 
q>0.3 r.l.u. As will be discussed in a later section, the lo- 
cal frequency minimum in the b-polarizcd branch is seen 
only at low temperatures. The dispersion of the bond- 
stretching modes with Ai-symmetry is, however, some- 
what obscured by an anticrossing with a c-axis polarized 
branch which starts from the Ag mode at 62 meV. In 
the b* direction, a further complication arises from the 
presence of another branch of the same symmetry in the 
energy region of interest: the corresponding zone-center 
mode has an energy of E=59 meV and involves longitu- 
dinal Cu-0 chain bond-stretching vibrations. When ap- 
proaching the zone boundary, there is another problem 
for the determination of the bond-stretching phonon fre- 
quencies: the intensities associated with bond-stretching 
modes overlap with those from other branches dispers- 
ing upwards in the energy range 45-55 meV (Figs. 4,5). 
They have, in principle, in-plane oxygen bond-bending 
character, but there will be inevitably some hybridiza- 
tion with the bond-stretching mode once the two modes 
are very close in energy. We note that this hybridization 
vanishes only close to the zone boundary, because there 
the bond-bending modes and the bond-stretching modes 
belong to different symmetry classes. 

Motivated by reports [15j of an unusual dispersion of 
transverse bond-stretching branches in YBa2Cu307 we 
studied these branches as well. As can be seen in Fig. 
9, our data do not support the claim that the transverse 
branch along b* has a rather strong dispersion up to 
q=0.25; to the contrary, this branch has as little disper- 



sion as the corresponding one along a* 20] . Figure 9 also 
shows the TO branch starting from the longitudinal chain 
O mode. There was no problem to follow this branch up 
to the zone boundary but attempts to follow the corre- 
sponding branch in the longitudinal direction had only 
very limited success. We learned from our model cal- 
culations that the LO branch starting from the chain O 
mode hybridizes very strongly with the other branches of 
the same symmetry - in particular, with the c polarized 
apical O vibrations - already at small wave vectors. For 
this reason, we were unable to follow this branch beyond 
q=0.15 r.l.u. 



B. Temperature dependence 

Selected phonons were studied as a function of temper- 
ature on a twinned sample. We found that the frequen- 
cies of the zone center phonons and of the zone bound- 
ary phonons show very little change when going from 
T=12 K to room temperature (Fig. 10). A certain in- 
crease in linewidth as is evident for the zone boundary 
phonon is typical of anharmonic behavior. The decrease 
of the phonon intensities with increasing temperature can 
be accounted for by the Debye- Waller- factor. Our find- 
ing that even the strongly renormalized zone boundary 
phonons show very little temperature dependence was in 
conflict with claims made by Chung et al.^f| that the 
bond-stretching phonons show an anomalous tempera- 
ture behavior throughout the Brillouin zone (see Fig. 9 
of Ref. 15). Therefore, we finally embarked on a de- 
tailed study to check these claims. Although we could 
not confirm many of the temperature induced changes 
reported in Ref. 15 we found indeed a very pronounced 
effect of certain phonon modes. As explained in a sepa- 
rate publication 17] a downward shift of spectral weight 
by at least 10 meV was observed for b polarized bond- 
stretching phonons of A4 '-symmetry at qw0.27 r.l.u. on 
cooling from room temperature to T=12 K. The tem- 
perature evolution starts well above the superconducting 
transition temperature, i.e. at about T=200 K. Data for 
T=200 K have been included in Fig. 1. These data were 
important for adjusting the parameters of our model. 
When the parameters were adjusted to better fit the data 
of the A4 '-branches, we obtained an improved fit for the 
Ai '-branches as well, which further validates our analy- 
sis. 



V. DISCUSSION 

The present results qualitatively confirm the results of 
an early study by Reichardt^4| but also go much beyond 
it in that we were able to elucidate the a-b anisotropy 
of the bond-stretching modes. Furthermore, we have 
gained some insight into the origin of the broad inten- 
sity distributions observed in the region half-way to the 
zone boundary. It appears that that the complex line- 
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shapes observed for a-polarized bond-stretching modes 
of Ai-symmetry can be attributed in large part to anti- 
crossings with other branches of the same symmetry. We 
agree with the observation reported in Ref. 15 that the 
bond-stretching branch along a* splits into two branches 
half-way to the zone boundary; however, we disagree on 
the interpretation. The phenomenon is the natural con- 
sequence of hybridization between the O in-plane bond- 
stretching modes and those of apical O vibrations along 
c. This phenomenon is illustrated in Fig. 11 for the zone 
boundary point. As for the b direction, there are two rea- 
sons for the observed complex line shapes: (i) anticross- 
ings of different branches and (ii) the anomalously steep 
dispersion - probably associated with very broad intrinsic 
linewidths - at q half-way to the zone boundary. In order 
to study the importance of (ii) we simulated the scans 
shown in Figs. 3,5 using the model fitted to the 200 K 
data, i.e. a model giving a smooth dispersion throughout 
the Brillouin zone. The following conclusions are drawn 
from this simulation (Fig. 12) Although the differences 
between Fig. 12 and Fig. 5b are somewhat subtle, there 
is definitely a better agreement between the experimen- 
tal results displayed in Fig. 5a and the calculated ones in 
Fig. 5b than with those in Fig. 12. This indicates that 
the bond-stretching branch of Ai '-symmetry exhibits the 
same type of anomalous low-T behaviour as its counter- 
part of A4 '-symmetry. This conclusion is in line with 
evidence for an anomalous temperature dependence of 
A 1 '-modes reported in Ref. 15. On the other hand, it 
would have been extremely difficult to establish details of 
the temperature dependence from measurements of the 
A 1 '-branches alone in view of the formidable complica- 
tion due to the anticrossings. Twinning of the sample 
further aggravates the situation considerably. Therefore, 
it is not surprising that the temperature dependence re- 
ported in Ref. 15 agrees with our results only on a qual- 
itative level. 

Measurements of the O in-plane bond-stretching 
phonons in optimally doped Lai.ssSro.isCuC^ revealed 
a considerable broadening of the phonon lines for q>0.1 
r.l.u..Q, 0. This broadening was found to be partic- 
ularly pronounced around q=(0.3,0.0) 2j. We expect a 
similar effect for optimally doped YBCO. Unfortunately, 
the problems for extracting linewidths discussed above 
do not allow us to deduce intrinsic phonon linewidths 
from the available data. We hope to achieve progress on 
this matter by a study of the A4 resp. the A4' modes on 
a detwinned sample. 

Our model calculations indicate that the observed a- 
b anisotropy of the bond-stretching branches can be 
largely attributed to the orthorhombicity of the struc- 
ture (apart from the anomalous temperature dependence 
along b*)[22|]. That is to say, the general energy differ- 
ence between the two directions is reproduced by a simple 
potential model in conjunction with the experimental dif- 
ference between the a and b lattice parameters. Further, 
the special terms needed to lower the zone boundary fre- 
quencies are of similar size. This contrasts with the re- 



sults of a previous study of underdoped YBa 2 Cu306.6 
which had shown a much stronger a-b anisotropy at the 
zone boundarv|27j which remains to be understood. 

We have already mentioned in Sec. Ill the results of a 
recent ab-initio calculation |23( of the phonon dispersion 
in YBa2Cu307 using Density Functional Theory in the 
local density approximation. The theory quantitatively 
predicts the frequency drop from the zone center to the 
zone boundary for the bond-stretching modes both along 
a* and b* at T=200 K; however, it also predicts a con- 
siderable downward dispersion for the [110] direction as 
well («8 meV) which is not borne out by experiment [l4|. 
(Also, it is not possible to check the doping dependence 
because of the known limitations of LDA for describing 
the correlated-insulator parent compound, YBa2Cu306-) 
Turning to a different mode, the theory makes an accu- 
rate prediction of a rather low frequency for the longi- 
tudinal chain O vibrations (61.5 meV). This is not quite 
as low as observed in experiment |28j (59 meV) but much 
lower than expected from the simple potential modelfl8| 
(75 meV) in which the Cu-0 potential was assumed to be 
the same for the Cu-0 chains and for the Cu-0 planes. 
Unfortunately, we were unable to check the theoretical 
predictions that the longitudinal branch starting from 
the chain O mode has a strong downward dispersion and 
exhibits a local minimum at about half-way to the zone 
boundary, related to Fermi surface nesting of chain O 
electronic states. The strong mixing of the chain O mode 
with other modes of the same symmetry did not allow us 
to isolate the contributions of the chain O vibrations from 
the rest. 

We note that in another cuprate, i.e. (La,Sr)2Cu04, 
the phonon renormalization upon doping is also much 
stronger in the [100] direction than in the [110] 
direction^]] . Thus, this behavior seems to be common 
to Cu-0 planes, but independent of extraneous features 
such as chains. The doping dependence and spatial 
anisotropy of the phonon softening appear to be captured 
by a calculation in which the t-J model is extended to ex- 
plicitly include electron-phonon couplings. 0, Ell- R6sch 
and Gunnarsson[24| have evaluated the doping and q de- 
pendence of the O bond-stretching mode of a Cu-0 layer 
through a t-J model with electron-phonon interactions 
derived from a three-band model. With their purely two- 
dimensional model, they obtain the correct anisotropy 
between the [1,1] and [1,0] directions, and reasonable re- 
sults for the doping- and q-dependent softening of the 
mode along the [1,0] direction. 

From the discussion above, it is evident that the lat- 
tice dynamics of YBa2Cu307 is less anomalous than ad- 
vocated in Ref. 15, but nevertheless we think that the 
behavior of the bond-stretching modes indicates a strong 
electron-phonon coupling. The question of whether or 
not this electron-phonon coupling is relevant for high 
T c superconductivity cannot be directly answered from 
the phonon data alone. From their theoretical analysis, 
Bohnen, Heid, and Krause[2^ came to the conclusion 
that conventional electron-phonon coupling is much too 
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weak to explain a high T c , in contrast to the case of 
MgB 2 where conventional electron-phonon does seem to 
explain the observed T„|2flj|. 



VI. CONCLUSIONS 

The present study has elucidated the a-b anisotropy 
of the plane-polarized bond-stretching vibrations in op- 
timally doped YBa 2 Cu3 06.95. We confirmed that the 
bond-stretching branches show a steep downward disper- 
sion in both the a and the b directions at q«0.25 r.l.u. 
that leads to anti-crossings with c polarized branches. 
The resulting anti-crossing gaps can be explained by lat- 
tice dynamical calculations using a shell model. In the a 
direction, the bond-stretching mode of Ai-symmetry is 
strongly hybridized with a c-polarized branch over a large 
part of the Brillouin zone, which explains the "splitting" 
of this branch reported in Ref. 15 in a natural way. The 
dispersion of the bond-stretching branches is very well re- 
produced by density-functional theoryj^ if the compari- 
son is based on 200 K data, whereas the pronounced soft- 
ening observed in b-polarized branches at low T remains 
unexplained. Although the pronounced renormalization 
of zone boundary bond-stretching modes on going from 
insulating 06 to superconducting 07 indicates a strong 
electron-phonon coupling, theory 23] suggests that it is 
too weak to be relevant for high T c superconductivity. 
Thus, the temperature effect, presented in detail in Ref. 
17, seems to be the most anomalous phenomenon in the 
phonon properties of optimally doped YBCO. 
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APPENDIX A: APPENDIX A: ANTI-CROSSING 
BEHAVIOR 

YBa2Cu307 has 13 atoms in the unit cell and hence 
there are 39 branches of the phonon dispersion in each 
direction. These branches can be grouped according 



to symmetry of the eigenvectors of the atomic displace- 
ments. For the case of the [100] or the [010] direction, 
there are four different symmetry classes. Within each 
subset of dispersion curves, the polarization may be very 
different. Still, none of the branches is allowed to cross 
another one of the same subset. If, as a result of the spe- 
cial nature of the interatomic interactions, two branches 
tend to cross each other in a particular part of the Bril- 
louin zone, this will lead to a mixing of the polariza- 
tion patterns and to a mutual repulsion of the eigenfrc- 
quencies (exactly as in the problem of coupled harmonic 
oscillators). However, the degree of hybridization and 
the strength of the repulsion cannot be predicted from 
general principles but depend on the symmetry of the 
crystal lattice and the peculiarities of the interatomic 
interactions. For illustrative purposes, we consider the 
following case: a simple interatomic potential modelfisl 
predicts a flat dispersion for the two topmost branches 
of Ai-symmetry in the a direction (Fig. 13). These 
two branches have in-plane polarization and c-axis po- 
larization, respectively. As a consequence, the inelas- 
tic scattering structure factors calculated for momentum 
transfers along the line (3+x,0,0), i.e. in the basal plane, 
will be very different for the two branches. In a corre- 
sponding neutron scattering experiment, only the upper 
branch will be detected. The situation dealt with so far 
corresponds to the situation found in 06. In 07, how- 
ever, the upper branch acquires strong dispersion from 
electron-phonon coupling effects as discussed in this pa- 
per. For the sake of simplicity, we have simulated such 
effects by adding a single special term to the dynamical 
matrix lowering the frequency of only the bond-stretching 
vibrations. The term is designed to produce a maximum 
effect at q=0.25 r.l.u. (for details, see APPENDIX B). 
This term was tuned to reduce the frequency of the bond- 
stretching vibrations at q=0.25 to below that of the api- 
cal O vibrations. As a consequence, phonons of the lower 
branch acquire some a-axis polarization over an extended 
range of wave vectors making them observable at momen- 
tum transfer along a (Figure 13). Evidently, hybridiza- 
tion and repulsion of the two types of vibrations are quite 
strong in YBa2Cu307. 



APPENDIX B: APPENDIX B: SPECIAL FORCE 
CONSTANTS 

As mentioned in Section III. 1 , several special terms 
have to be added to the shell model to arrive at a quan- 
titative description of the bond-stetching branches. The 
term describing screening by free carriers has been ex- 
plained in detail in Ref. 30. The term accounting for 
a planar breathing deformability has been given in Ref. 
30 as well. The parameter A for this planar breathing 
deformability was tuned to reproduce the observed |L4j 
bond-stretching phonon dispersion in the [110] direction. 
Although the bond-stretching modes in the [110] direc- 
tion do show a considerable renormalization with the 
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metal-insulator transition, the branch remains flat. As 
a consequence, the parameter A was found to be rather 
small (A=25000 dyn/cm). 

The further terms needed to arrive at a satisfactory 
description of the bond-stretching branches in the [100] 
and in the [010] directions are ad hoc extensions of the 
term for the breathing deformability: a 2x2 sub-matrix 
is added to the dynamical matrix for the plane oxygen 
atoms 02 and 03 (and, of course, for the other plane oxy- 
gen atoms in the bi-layer as well) . The diagonal elements 
of the matrix are given by 

D 2 (l, 1) = B wo x sin 2 (nq x )+C W o-sin 2 (2nq x ))-\sin 2 (irqa;) 



D 2 (2,2) = B io* sin 2 (nq y )+C 0W -sin 2 (2-n-q x ))-\sin 2 (Tvq y ) 

The non-diagonal elements of the matrix are zero. 
For negative values of Bioo the frequencies of the bond- 



stretching modes are lowered in the [100] direction in a 
sine-like manner from the zone center to the zone bound- 
ary. On the other hand, this term has no effect in both 
the [010] and in the [110] directions. From a fit to the ex- 
perimental data we obtained Bioo=-70000 dyn/cm. The 
parameter Cioo describes the second harmonic of the lin- 
ear breathing deformability with a maximum effect at 
q=0.25 a* and zero effect at q=0 and at q=0.5 a*. The 
fit resulted in a value Cioo=-60000 dyn/cm. The corre- 
sponding terms for the [010] direction were obtained as 
B 010 =-75000 dyn/cm and C O i O =-55000 dyn/cm. 

In order to describe the anomalous softening upon 

■ 2/ Ml/2 

sl %3mS)g\ in the [010] direction at least in a semi- 
quantitative manner another term was added to D2(2,2) 
in analogy to the term with the prefactor Coio but with 
. the factor^in (2ir qy) raised to the nth power. This term 

81 waitlisted by trial and error to C' O io=-100000 dyn/cm 
and n=6. 
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FIG. 1: Illustration of the effect of the special terms in- 
cluded in the model calculations. Black line: no extra terms. 
Green line: after inclusion of a (negative) planar breathing 
deformability chosen such as to reproduce the experimen- 
tal zone boundary frequency in the (110) direction. Violet 
line: after inclusion of a special term to lower the frequency 
of the half-breathing' mode. Red line: after inclusion of a 
special term to lower the frequency of the bond-stretching 
modes halfway to the zone boundary. Blue line: after inclu- 
sion of a further term to simulate the anomalous softening 
at low temperatures. Only the branch with in-plane bond- 
stretching character is shown for clarity. Note that the special 
terms directly affect only the vibrations with bond-stretching 
character. Blue/Red circles denote phonon peak positions at 
T=12K/200K. The blue bar at q=0.275 r.l.u. denotes an ill- 
defined phonon peak containing also contributions from the 
next lower phonon branch with bond-bending character. 



FIG. 2: Measured (a) and calculated (b) neutron scattering 
intensity along the a direction from Q=(4,l,0) to (3.5,1,0). 
This Brillouin zone maximizes the scattering cross section 
for phonons with strong longitudinal plane oxygen vibrations 
perpendicular to the chains. Further, it allows one to ex- 
ploit focusing effects for branches with a downward disper- 
sion. Background, linear in energy and independent of Q, 
has been subtracted from the experimental raw data to em- 
phasize one-phonon scattering. The calculated one-phonon 
neutron scattering intensity was corrected for spectrometer 
resolution. 



FIG. 3: The same as Fig.l but for the b direction. The mea- 
surements were carried out from Q=(l,4,0) to (1,3.5,0). 



FIG. 4: Coulour coded contour plots of the measured and 
calculated neutron scattering intensities shown in Fig. 1. 
No data is available for the regions colored in brown. The 
calculated phonon dispersion curves (without the resolution 
correction) are shown as blue lines. According to the model, 
phonon branch anticrossing causes dramatic changes in eigen- 
vectors within the same branches between the zone center and 
the zone boundary, which results in corresponding changes in 
one-phonon scattering cross sections. Apart from this effect, 
there is a 30 % increase of the structure factors of the bond- 
stretching modes from the zone center to the zone boundary 
in the (4,1,0) Brillouin zone. 



FIG. 5: The same as Fig. 3 but for the b direction. The 
measurements were carried out from Q=(l,4,0) to (1,3.5,0). 



FIG. 6: Calculated dispersion curves of phonon branches of 
Ai-symmetry and of A 1 '-symmetry, respectively, as well as 
experimentally measured phonon peak energies. Arrows point 
to Ag modes, the other modes being of (B2 U ) character. 
Data in the energy range 10 meV<E<42 meV were taken 
from previous measurements on twinned samplesQ]. Data for 
the longitudinal acoustic branches were taken from previous 
measurements on a small (12 mm 3 ) detwinned sample 19]. 



FIG. 7: Neutron spectra taken on a twinned sample at T=12 
K with a high resolution configuration. The spectra were fit- 
ted with two resp. three Gaussians and a sloping background. 
The peaks are assigned to branches of A4-symmetry in the a 
direction or the b direction, respectively. 



FIG. 8: Dispersion of high-energy phonon branches of A4- 
symmetry along the a direction as determined on the twinned 
sample. The full lines were calculated from the model de- 
scribed in the text and the dashed lines were taken from Ref. 
23. Note that the ab-initio resultsfzj^ were shifted upwards 
by 1 meV. 



FIG. 9: Dispersion of the transverse bond-stretching branches 
along a* and along b*, respectively, as observed on an un- 
twinned sample at T=12 K. The high-energy branch and 
the low energy branch along a* have in-plane oxygen and 
chain oxygen character, respectively. Lines are a guide to 
the eye. The arrows denote the frequencies observed by 
ellipsometrv|2l| . 



FIG. 10: Energy scans taken on a twinned sample at two 
different temperatures. 



FIG. 11: Displacement patterns of the two highest energy 
zone boundary modes with wave vector q=(0.5,0,0) as calcu- 
lated from two different models A and B. Model A is described 
in the text and reproduces the anomalously low frequency of 
the in-plane Cu-0 bond-stretching mode by including special 
terms designed for that effect. This choice leads to a strong 
hybridization of the in-plane bond-stretching mode and of the 
c-axis polarized apical O mode with energies of 58 meV and 66 
meV, respectively. It corresponds to the situation displayed 
in Fig. 5. The special terms are ommitted in model B which 
leads to a much higher bond-stretching mode frequency (74 
meV) and consequently to a much smaller hybridization with 
the apical O mode (61 meV). 



FIG. 12: Colour coded contour plot of calculated 200K 
neutron scattering intensity for the line from Q=( 1,4,0) to 
(1,3.5,0). The calculated phonon dispersion curves along the 
high symmetry line are shown as blue lines. The model used 
for this plot differs from the one used for Fig. 5 by omit- 
ting the special term introduced to describe the anomalous 
softening observed at low temperatures. The corresponding 
behavior of the A4 '-branches is shown in Fig. 1 by a red line 
and a blue line, respectively. 



FIG. 13: Top: Calculated dispersion of the two highest Ai- 
branches before (full lines) and after (dashed lines) inclu- 
sion of a special term to reduce the frequencies of the bond- 
stretching modes half way to the zone boundary. Bottom: 
inelastic structure factors calculated along the line (3,0,0) to 
(3.5,0,0). 



TABLE I: Table I: Parameters of the shell model for 
YBa2Cu3C>7. Z(k), ionic charge; Y(k), shell charge; K(k), 
core shell force constant. The repulsive Born-Mayer potential 
used in Ref. 19 was replaced by force constants: F fcfc / and G^y 
are longitudinal and transverse force constants, respectively. 
Following Ref. 19, an attractive van der Waals potential with 
Cfcfc'=100 eV Awas assumed to act between oxygen atoms. 
Note that the chain copper and the chain oxygen are labelled 
as Cu(l) and 0(4), respectively. 
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